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A LOW-COST INERTIAL SMOOTHING SYSTEM 


FOR LANDING APPROACH GUIDANCE 

By Frank R. Niessen 
Langley Research Center 

SUMMARY 

Accurate position and velocity information with low noise content is required for 
instrument approaches and landings, for both control and display applications. However, 
current landing guidance systems provide position reference signals with so much noise 
that velocity information often cannot be obtained through differentiation without resulting 
in unacceptable noise or unacceptable time lags if filtering is employed to reduce the 
noise. An inertial navigation system or, perhaps, a Doppler velocity system could be 
used for this purpose, but such systems are expensive, particularly if redundant systems 
are required for a fail-operational capability. 

In a current VTOL automatic instrument approach and landing research program, 
landing guidance position reference signals have been smoothed by using a mix -filter 
technique. The inputs to the mix filter include both inertial acceleration inputs, which 
provide high-frequency position and velocity information, and also radar position inputs, 
which provide low-frequency position and velocity information. Low-cost aircraft instru- 
mentation, including attitude reference gyros and body-mounted accelerometers, have 
been used to provide the inertial acceleration information. Since the onboard inertial 
information is relied on for only short periods of time in this application, high-quality 
inertial sensors are not required. An in-flight comparison of signal quality and accuracy 
has shown good agreement between the low-cost inertial smoothing system and an aided 
inertial navigation system. Furthermore, the low-cost inertial smoothing system has 
been proven to be satisfactory in control and display system applications for both auto- 
matic and pilot -in -the -loop instrument approaches and landings. 

INTRODUCTION 

For the precision instrument approach and landing task, aircraft require accurate 
inertial (that is, ground referenced) position and velocity information. This information 
will be necessary for automatic landings to touchdowns and for time -constrained landing 
approaches in the terminal area. Inertial position and velocity information are particu- 
larly necessary for the V/STOL instrument approach and landing, because at low speeds 



these aircraft and, especially, their flight paths can be greatly affected by gusts and 
wind shears. It can be assumed that position information will be provided by a ground- 
based tracking station, such as the proposed microwave landing system, but this position 
information is expected to be too noisy to differentiate to obtain acceptable velocity infor- 
mation. An inertial navigation system or, perhaps, a Doppler velocity system could be 
used for this purpose, but such systems are expensive and if redundant systems were 
required, the cost could be prohibitive for most potential operators. 

A promising concept is that of a low-cost inertial smoothing system, where rela- 
tively low-cost aircraft instrumentation, attitude reference gyros and body-mounted 
accelerometers, can be used to provide short-term inertial navigation. This informa- 
tion is mixed with position information from the landing guidance system to obtain accu- 
rate estimates of velocity and position. This report presents the theoretical background 
which has been used in dealing with this particular estimation problem. In addition, 
flight-test results are presented for a low-cost inertial smoothing system which has been 
developed for use in a current VTOL automatic instrument approach and landing research 
program. 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements and cal- 
culations were made in U.S. Customary Units. 

A,B,C matrices 

a x ,a y ,a z body -mounted accelerometer outputs, m/sec 2 (ft/sec 2 ) 
b,c single column matrices or vectors 

g gravity constant, 9.8 m/sec 2 (32.2 ft/sec 2 ) 

K gain matrix 

k element of gain matrix 

R Riccati equation matrix solution 

s Laplacian operator 

t time, sec 
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u control input vector 

v input noise vector 

w measurement noise vector 

X,Y,Z displacement in runway reference coordinate frame (see fig. 3), m (ft) 

inertial accelerations in aircraft reference coordinate frame 
(see fig. 3), m/sec2 (ft/sec 2 ) 

x state vector 

y output vector 

£ damping ratio 

0 pitch attitude, positive nose upward, rad 

a y 2 variance of input signal noise 

cr w 2 variance of measurement signal noise 

r time constant, sec 

< p roll attitude, positive right wing down, rad 

yaw attitude, positive nose right, rad 
co n undamped natural frequency, rad/sec 

Superscript: 

T matrix transpose 

A dot over a symbol indicates a derivative with respect to time. A circumflex (*) 
denotes an estimated value. 
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THEORETICAL BACKGROUND 
Differentiation of Position Information 

There are inherent problems in differentiating position information to obtain veloc- 
ity information, because of signal noise. Of course, position and velocity information 
which is to be used in display or control system applications must be practically noise 
free. Figure 1 shows the transfer function for approximate differentiation, which is a 
typical approach. This method is referred to as approximate differentiation because the 
output is actually the rate of change of the input filtered with a first-order lag with a 
time constant r. The time response shows that for a unit step input, the resultant out- 
put is a pulse with amplitude 1/r. The effect of filtering or lagging a feedback signal is 
to reduce the stability of the closed-loop system. From this standpoint, it is necessary 
that t be as small as possible. However, it can be seen that any high-frequency noise 
which is part of the position signal will be amplified by a factor of 1 /t; thus, to attenuate 
the high-frequency noise, it is necessary that t be larger than 1.0. 

A trade-off has to be made between lag and noise level, but more frequently than 
not, an acceptable trade-off cannot be made. This is true simply because the noise level 
associated with position information from current landing guidance systems is such that 
even with the maximum acceptable lag, an acceptable level of noise for the velocity out- 
put cannot be obtained. For example, even if a lag of 0.5 second and a peak-to-peak 
noise level of 0.3048 m/sec (1.0 ft/sec) would be acceptable, these conditions would 
require a position signal with a peak-to-peak noise level of only 0.1524 m (0.5 ft), which 
is unrealistically low for current landing guidance systems. Figure 2 is a time history 
obtained from the GSN-5 precision radar facility at Wallops Station which is described 
in appendix A. The velocity information was obtained by the GSN-5 computer by use of 
an approximate differentiation circuit with a time constant of 0.5 second. The peak-to- 
peak noise on the velocity signal is on the order of 3.048 to 6.096 m/sec (10 to 20 ft/sec) 
and thus indicates that the peak-to-peak noise on the position signal is 1.524 to 3.048 m 
(5 to 10 ft). So, even for a precision tracking radar, the position signal noise is an order 
of magnitude greater than the value that would be acceptable. The proposed microwave 
landing system, as noted in reference 1, has a specification for a one -sigma accuracy of 
only 3.048 m (10 ft) for slant range for the most advanced configuration. For these rea- 
sons, differentiation of position information, by itself, does not appear to be a very 
promising approach. 


Inertial Smoothing Concept 

V 

The inertial smoothing mix filter combines acceleration with position data to deter- 
mine low-noise estimates of both velocity and position. The coordinate frames of refer- 
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ence which were used in measuring aircraft acceleration and position are depicted in 
figure 3. A simplified approach to the estimation problem was adopted by using three 
independent mix filters, any correlation in the noise characteristics of the X, Y, and 
Z position signals being ignored. Figure 4 illustrates the mix -filter configuration which 
was used for inertial smoothing. It is noted that the form of the mix filter for this par- 
ticular estimation problem is identical to that of a Kalman filter. The part of the system 
drawn with solid lines (fig. 4) represents the high-frequency computation of velocity and 
position, based on aircraft acceleration. There is no lag as acceleration is integrated 
directly once to obtain velocity and twice to obtain position. Note that the mix filter pro- . 
vides a position estimate as well as a velocity estimate. As drawn with dashed lines 
(fig. 4), the difference between the estimated position and the position measured by the 
tracking station is fed back as a correction to both the velocity estimate and the accelera- 
tion input. In contrast to approximate differentiation, the velocity estimate is the output 
of an integrator which attenuates noise on the position input and also the accelerometer. 
For additional information, table I contains transfer functions which indicate the response 
of each of the estimator outputs to individual acceleration and position inputs. 

The selection of the mix -filter gains was based on the general steady-state Kalman 
filter solution for this particular estimation problem. Given the noise properties of the 
inputs, the Kalman filter solution provides the gains for an optimal estimator in the sense 
that the estimates will have minimum variance noise. The results of the general solu- 
tion for this particular filter are discussed here, and a more detailed treatment of the 
solution is given in appendix B. 

The gains kj and k 2 can be expressed in terms of the familiar second-order 
parameters £ and w n as kj = 2£w n and k 2 = u> n 2. The damping ratio was found 
to be constant and the undamped natural frequency was determined to be a function only 
of the ratio of accelerometer noise to position noise 



w n - ^°v/ a w 

It is important to note that the solution is dependent only on the relative noise between the 
inputs and not on the absolute noise levels. It is possible that if the absolute noise levels 
were too high, the mix filter would provide estimates which could be unacceptably noisy. 

FLIGHT-TEST VALIDATION OF INERTIAL SMOOTHING CONCEPT 

System Description 

A system based on the concept of inertial smoothing has been used in a current 
VTOL automatic approach and landing research program. This system, which will be 
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described here, was installed in the research helicopter shown in figure 5. The onboard 
sensors used were low-cost autopilot quality sensors as compared with more expensive, 
higher quality sensors that are used in inertial navigation systems. Several approxima- 
tions were made to minimize computational requirements, based on the assumption of 
small pitch and roll angles during the final approach and landing maneuver for the 
research tests. These approximations did not degrade the system for its intended use, 
nor do they in any way imply a limitation of the concept itself. 

The general arrangement of the system is depicted in figure 6. The body-mounted 
longitudinal and lateral accelerometers have been corrected for the effects of gravity by 
using the sine of pitch and roll angles from the onboard vertical gyro 

Xh = ^ - g sin 8 

Y h = ay + g sin 0 

The normal accelerometer has been corrected for effects of gravity by assuming small 
pitch and roll angles 

Z h = a z - g cos 8 cos <f> 

Z h * a z - g 

In resolving the accelerations along the body axes into the runway reference coordinate 
frame, it has been assumed that the longitudinal and lateral axes are nearly in the hori- 
zontal plane so that 

X * X h cos {\J/ - <// 0 ) - Y h sin(i// - * 0 ) 

Y ~ Y h sin(i// - + Y h cos(^ - x // Q ) 

Z « Z h 

The horizontal accelerations were resolved by using a sine -cosine resolver driven by a 
directional gyro synchro output. A differential synchro input was incorporated to permit 
selection of any desired runway -reference heading The three position signals X, 

Y, and Z were obtained from the GSN-5 radar station by means of narrow-band fre- 
quency modulated (FM) telemetry equipment. The GSN-5 system itself is described in 
appendix A. A 10-volt analog computer has been used to perform all the necessary com- 
putations onboard the aircraft. Figure 7 is the analog computer schematic which shows 
the’ details of these computations and indicates the scaling which was used. 

If the accelerometers were slaved to the vertical so that the longitudinal and lateral 
accelerometers would indicate true horizontal accelerations and the normal accelerom- 
eter would indicate true vertical accelerations, then the need to correct the longitudinal 
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and lateral accelerometers for effects of gravity would be eliminated and would make the 
rest of the computations valid for other than small pitch and roll angles. This was not 
done, however, for the system described herein. 

The mix -filter gains which were used corresponded to a mix -filter natural fre- 
quency of co n = 0.45 rad/sec with a settling time constant of 12.5 seconds, based on the 
time to settle within 2 percent of steady state. These gains were selected on the basis of 
the reasoning that the time constant should be long enough to insure that noise from the 
radar position signal would be satisfactorily attenuated, but short enough to insure that 
errors which would result from inaccuracies associated with the acceleration information 
would be kept small. 


Flight-Test Results 

As noted previously, the inertial smoothing system was developed for use in a cur- 
rent VTOL instrument automatic approach and landing research program. The modifica- 
tion of the recording system that was necessary to obtain the data presented here was 
restricted in order that the data could be obtained in a timely manner, without impeding 
the main research program. For this reason, data were obtained for only X and Y in 
one instance, and for only X in another instance. Nevertheless, the axes that were 
selected for documentation were the axes with the least desirable scale factors and, con- 
sequently, represent the worst case rather than the best. 

Figures 8 and 9 show the input and output signal noise characteristics for the iner- 
tial smoothing system. These data were recorded separately for X and Y during 
hovering flight near the landing zone. These data were obtained by an FM magnetic tape 
recording system and were later sampled at a rate of 100 Hz during the data -reduction 
process. Longitudinal acceleration was not recorded during these tests; however, its 
signal noise characteristics are very similar to those of the lateral accelerometer. The 
noise on the position output as compared with the noise on the position input signal, which 
has a peak-to-peak amplitude of approximately 30.48 m (100 ft), has been greatly reduced. 
The noise on the position input signal was mostly due to telemetry noise, which was nearly 
1 percent of full scale. The accelerometer signal noise, mainly due to aircraft structural 
vibration, has also been essentially eliminated by the integration process, as evidenced 
by the velocity estimate. It was determined that the signal noise contributed by the analog 
computer components themselves, which depends on computer scaling, would be expected 
to result in approximately 0.061 m/sec (0.2 ft/sec) peak-to-peak noise for the velocity 
estimate for X, for which the scaling problem was most severe. This level of computer- 
generated noise accounts for nearly all the noise observed on the velocity output; hence, 
the mix filter has essentially eliminated the effects of both acceleration and position mea- 
surement noise. 
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A comparison in accuracy was made in flight between the inertial smoothing system 
and an aided inertial navigation system. The aided inertial navigation system was a modi- 
fied Gemini inertial platform system. Baseline performance data on the modified system 
are presented in reference 2, whereas details of this system’s navigation computations, 
including the update logic, are described in reference 3. Both systems relied on the 
GSN-5 radar for long-term position and velocity information; the inertial smoothing sys- 
tem received continuous position information, whereas the aided inertial navigation sys- 
tem received position updates at 1.0-second intervals. Without updates, the position- 
error drift rate of the inertial navigation system was approximately 2.0 nautical miles 
per hour, whereas the position-error drift rate of the inertial smoothing system, without 
position feedback, was estimated to be on the order of 100 to 200 nautical miles per hour. 
This high drift rate was mainly due to approximations which were made in resolving the 
accelerations and also was a result of computer scaling limitations. Figure 10 shows 
close agreement between the outputs of the two systems for a decelerating approach to 
hover. Figure 11 is a hovering translation time history which compares the two systems 
on a much more sensitive scale. Agreement between the velocity outputs is within 
0.3048 m/sec (1.0 ft/sec). 

In addition to the tests discussed, the six position and velocity outputs of the inertial 
smoothing system have been used in both control and display applications in a VTOL auto- 
matic approach and landing research program, for which the system was developed. Dis- 
plays driven by these outputs consisted of a flight -director indicator, a horizontal- 
situation moving-map display, lateral and vertical flight -path error needles, and simu- 
lated radar altimeter. The guidance computer, which provided pitch, roll, and power 
flight director display commands, also provided similar commands to the control system 
in the automatic approach mode. The VTOL landing approach task in these tests, involved 
acquisition of the runway center line, capture of a 6°, 10°, or 15° glide path, deceleration 
to a hover, vertical descent, and touchdown. The entire sequence could be accomplished 
either manually (with the pilot centering the flight director commands) or automatically. 
Signal noise was not apparent either in the display movements or, while in the automatic 
mode, in the control actuator motions. The accuracy of the inertial smoothing system 
enabled tracking of the approach path and speed profiles with a very high degree of pre- 
cision. The manual and automatic VTOL approach results obtained by using the inertial 
smoothing system are described in reference 4. 

CONCLUSIONS 

A study has been conducted at the Langley Research Center to investigate the 
feasibility of using low-cost, conventional aircraft instrumentation in combination with 
landing guidance system signals to provide acceptable position and velocity information 
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for landing approach guidance. Based on the work described in this report, the following 
conclusions have been drawn: 

1. Differentiation of landing guidance system signals, by itself, will not provide 
adequate velocity information for use in either controls or displays since differentiation 
is a process which, by its very nature, tends to amplify noise. Current landing guidance 
systems do not provide signals which are sufficiently noise free to permit differentiation, 
nor does it appear likely that systems, which are now being proposed, will be capable of 
providing such noise -free signals. 

2. Flight-test data indicated that accelerometer and radar position signal noise, 
approximately 1.2192 m/sec 2 (4 ft/sec 2 ) and 30.48 m (100 ft) double -amplitude noise, 
respectively, had been greatly reduced by the inertial smoothing mix filter, which pro- 
vided a velocity estimate with only 0.06096 m/sec (0.2 ft/sec) double -amplitude noise. 

3. An in-flight comparison of signal quality and accuracy between the low-cost 
inertial smoothing system and an aided inertial navigation system showed agreement to 
within 0.3048 m/sec (1.0 ft/sec) for the velocity outputs. 

4. The concept of using low-cost autopilot quality aircraft instrumentation to pro- 
vide inertial smoothing of landing guidance system signals has been validated through the 
use of such a system in a VTOL automatic approach and landing research program, for 
which the system was developed. The VTOL landing approach task involved acquisition 
of the runway center line, capture of a 6°, 10°, or 15° glide path, deceleration to a hover, 
vertical descent, and touchdown; thus, the possibility of broad applications for this con- 
cept is indicated. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 18, 1973. 
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APPENDIX A 


GSN-5 PRECISION RADAR 

Position was sensed by the GSN-5 precision tracking radar system, located at 
Wallops Station, Virginia, where the flight tests were performed. A photograph of the 
GSN-5 is shown in figure 12. The position of the aircraft is sensed directly in terms of 
slant range, and azimuth and elevation angles of the radar antenna. This information 
is transformed into rectangular coordinates in the runway reference frame and trans- 
mitted to the aircraft by means of a narrow -band frequency modulated (FM) telemetry 
link. 

The radar is K-band and has an antenna beamwidth of approximately 0.5°. A pas- 
sive reflector has been mounted on the nose of the aircraft to provide a specific point of 
high-energy return in order to prevent skin tracking. The limits of the radar tracking 
antenna are 0° to 30° in elevation and ±45° in azimuth. The accuracy of the radar as 
specified by the manufacturer is 0.02° for the azimuth and elevation angles and 3.048 m 
(10 ft) or 1 percent (whichever is greater) for slant range. 
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APPENDIX B 


FIXED -GAIN SOLUTION FOR INERTIAL SMOOTHING MIX FILTER 

The general Kalman filter solution, obtained from reference 5, is outlined below. 

The plant dynamics are expressed as: 

x = Ax + Bu 

with measured outputs 

y - Cx + w 

The input and output measurement noise characteristics are specified by the correlation • 
matrices Q and P, respectively. The Kalman filter gain matrix, for the stationary 
case , is found from 

K = R 0 C T P -1 

where R 0 is the steady-state solution of the matrix Riccati equation 
R= AR + RA T - RC t P“ 1 CR + BQB t 

For the inertial smoothing mix -filter problem, as shown in figure 13, the accelera- 
tion input is regarded as a control input and, as a single input, is a scalar; therefore, 
the B matrix is reduced to a column matrix or vector b. Similarly, there is only one 
measured output; therefore, the C matrix becomes a row matrix c**. 

x = Ax + bu 

and b = 

T 

y = cx + w 

where c = 

The correlation matrices Q and P are reduced to the scalar quantities a 2 
and a w ^. The matrix Riccati equation, therefore, becomes 

R = AR + RA T - Rc c T R + b (oy 2 )b T 

• • 

The steady-state solution is found by setting R = 0. Taking each element of R, 

■ 0 - r21 + ri2 ■ {^ u2 




where A 


_[0 1 
~ _° °_ 
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APPENDIX B - Concluded 


ri2 = ° = r 


22 " ^f 11 ' 12 


r 2i = ° = r 


22 '1^2) r2iri1 


r 22 = 0 = ~( -Mr«, r, n + a. 




2/ 21 12 T u v 


By using the fact that the R matrix is symmetric and that r 2 j = r^, these equations 
can be reduced to 

/ 1 \_ 9 . 


0 = 2r 
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•'w 


2*11 


0 = r 


22 




~2 r 12 r ll 


0 

From these equations, the steady-state values for the elements of R have been found 
to be 


R 0 = 


r ll 

r 12 


^c v >/2 aw 3/2 

1 

r 21 

r 22 

11 

°v a w 

S„ 3 / 2 a 1/2 

» V w 


The gains are obtained from 



and result in kj = \J2 \ja v /a vf and k 2 = CT v/ a w' By ex P ressin g the gains as kj = 2t> n 
and k 2 = w n 2 , it may be shown that w n = \jo v /a w and £ = V 2/2 = 0.707. 
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TABLE I.- TRANSFER FUNCTION RELATIONSHIPS FOR 


INERTIAL SMOOTHING MIX FILTER 


Input 

Output for - 

• 

X 

X 

X 

s + kj 

1 

s 2 + k^ + kg 

s 2 + kjS + kg 

X 

k 2 s 

skj + kg 

s 2 + kjS + kg 

s 2 + kjS + kg 




X (approximate) 






+ 


t/i 


H 

X 
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Figure 4.- Inertial smoothing mix filter. 
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Figure 6.- General arrangement of inertial smoothing system. 
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Figure 8.- Inertial smoothing system input and output signal noise characteristics X. 

(Note that a x is not recorded.) 
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Figure 9.- Inertial smoothing system input and output signal noise characteristics Y. 

(Note that a x is not recorded.) 
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Figure 10.- Comparison of low-cost inertial smoothing system and an aided inertial 

navigation system. Decelerating approach. 
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Figure 11.- Comparison of low-cost inertial smoothing system and an aided inertial 

navigation system. Hovering flight. 







Figure 12.- GSN-5 precision-tracking radar facility. 
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Figure 13.- State -variable representation of inertial 
smoothing mix -filter problem. 
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